arXiv:math-ph/0112020v5 28 Feb 2003 


Transform of Riccati equation of constant coefficients through 

fractional procedure 

H.C. Rosu^*, A.L. Maduefio^, J. Socorro^ 

^ Dept, of Appl. Math., IPICyT, Apdo Postal 3-74 Tangamanga, San Luis Potosi, MEXICO 
^ Instituto de Fisica, Universidad de Guanajuato, Apdo Postal E-ljS, Leon, MEXICO 
We use a particular fractional generalization of the ordinary differential 
equations that we apply to the Riccati equation of constant coefficients. By 
this means the latter is transformed into a modified Riccati equation with the 
free term expressed as a power of the independent variable which is of the 
same order as the order of the applied fractional derivative. We provide the 
solutions of the modified equation and employ the results for the case of the 
cosmological Riccati equation of FRW barotropic cosmologies that has been 
recently introduced by Faraoni. 
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I. INTRODUCTION 

The fractional calculus is a generalization of the ordinary differential and integral calculus 
[1]. The main point is how to think about the derivative of order r, where r is an arbitrary 
real or complex number. In 1695, L'Hospital was the first to ask in a letter to Leibnitz on 
the possibility to perform calculations by means of a fractional derivative of order r = ^. 
Leibnitz answered that the question looked as a paradox to him but he predicted that in the 
future useful consequences might occur. In 1697, Leibnitz referring to the inhnite product 
of Wallis for 7r/2 used the notation and surmised that the fractional calculus could be 
used to get the same result. 

In 1819 the hrst mention of derivatives of arbitrary order occurred in a published text. 
The French mathematician S.F. Lacroix published 700 pages on the differential calculus, 
where one can find less than 2 pages dedicated to the fractional topic, which seems to be 
based on a result of Euler dated 1730. He started with 
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= X 


( 1 ) 


where n is an integer and wrote the mth derivative in the form 
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Next, Lacroix changed the factorial using the F(a:) function (introduced by Legendre), n 
was changed from an integer to a real number, n ^ a, and m was chosen to be m = A thus 
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^ r(a + l) 1 

dx^/"^ r(a + ^) 


( 3 ) 


In this way, he expressed the derivative of order one half by an arbitrary power of x. A 
simple example given by Lacroix refers to the case y = x 


ddl'^x 2^/x 


( 4 ) 


Along the years, great mathematicians, such as Euler, Fourier, Abel and others, did some 
work on the fractional calculus that nevertheless remained as a sort of curiosity. 

The modern epoch started in 1974, when a consistent formalism of the fractional calculus 
has been developed by K.B. Oldham and J. Spanier [2], It has been found that Lacroix’s 
result (4) coincides with that obtained by means of the present dehnition of the Riemann- 
Liouville fractional derivative. 


II. BASIC DEFINITIONS OF THE FRACTIONAL CALCULUS 


(i) One can dehne the fractional integral of order a > 0 as follows 

1 U f{x')dx' 


Xo 




r(Q;) JxQ {x — x'Y “ 

In particular, for Xq = 0 one usually writes 

’ r(a)Jo (x-x'Y-“' 

(ii) For /5 > 0 one can dehne the fractional derivative of order f3 in the following way 
dxd ^ dx"^ ^ dx'^T{n — j3) Jo {x — x'Y~'^~^d 


( 5 ) 


( 6 ) 


( 7 ) 


where n > j3. Thus the (3 fractional derivative is dehned as an ordinary derivative of order 
n of the fractional integral of order n — 13. 

(hi) The chain rule has the form 


^/(sW) = E (f), (^i) ;^/(9W). 

where k ^ N and are the coefficients of the generalized binomial 

r(i + /5) 


(0 


r r(l + A;)r(l - A; +/3) 


( 8 ) 


(9) 


(iv) Leibnitz’s rule for the derivative of the product has the form 


dxh 


oo Jk nP—k 

f{x)g{x) = YJ , 


k=0 
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where k E N. 


( 10 ) 
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III. ORDINARY DIFFERENTIAL EQUATIONS FROM THE FRACTIONAL 

CALCULUS 


Let L be a differential operator and let its action L[u{x)] = g{x) be used to dehne a 
differential equation. We generalize L to the fractional calculus by means of the fractional 
derivative of order = 1 — 5, 5 G (0,1), writing the following fractional differential equation 

Dl ^L[u{x)] = g{x). 

(11) 

When 5 = 1, then Dl~^L[u{x)] = L[u{x)]. 

One can obtain a solution of (11) by applying the fractional integral of the same order 
to the left 

{dI-%u{x)]} = D-^^-^^[g{x)], 

(12) 

L[u{x)] = D~^^-^'>[g{x)], 

(13) 

L[u{x)] = Di-^[g{x)]. 

(14) 

For example: L = -^ + p{x) implies 


L[u(x)] = -^u(x) + p(x)u(x), 
dx 

(15) 


■^u{x) +p{x)u{x) = Di ^[^(x)], (16) 

then 

+ c , (17) 

where 

jj,{x) = exp p{s)ds^ . (18) 

Equations (17) and (18) provide the solution to the linear generalized equation of the hrst 

order. This is not the unique possible fractional generalization [3,4], One could have taken 
Dlr^u{x) +p{x)u{x) = g{x) or some other procedure. However, the present approach leads 
to analytical results in applications. 
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IV. APPLICATION TO THE RICCATI EQUATION OF CONSTANT 

COEFFICIENTS 


The ordinary Riccati equation of constant coefficients is + av?‘{x) = 6, where a and 
b are constants. Thus the operator of the Riccati type is Ln = -^ + au{x) acting always in 
the space of functions u{x), i.e., Lru{x) = + av?‘{x) = g{x). Fractional considerations 

related to this operator can be found in the work of Metzler et al [4]. 

The fractional Riccati equation according to the scheme proposed in the previous section 
is 


D 


l-<5 

X 


A. 

dx 


+ au{x) 


u{x) = b, 


(19) 


or 

AAA. + au^x) = 

dx \ J X 

The right hand side can be written 

I|n-(1-S),l ^ '> r ^ ^ 

‘ ^ r{i-S)Jo (x-t)‘ r(i-^)(i-i) r(2-^)' 


( 20 ) 


( 21 ) 


Thus, solving the fractional Riccati equation is equivalent to solving the following type of 
particular, ordinary Riccati equation that we call the (5-modi£ed Riccati equation 


du{x) 

dx 


+ a[u{x)]‘^ 


bx^ ^ 


( 22 ) 


A. Solution of the (5-modified Riccati equation 


In order to solve (22) we use the transformation u = y'/ay^ u' = [y''y — {y'Y]/[ay‘^]. One 
gets the associated linear second order differential equation 


(Py{x) 

dx"^ 


ab 


x^ ^y = 0 . 


(23) 


Multiplying by x^ leads to 


X 


> d^y{x) 

dx"^ 


ab 


r(2 - h) 


x^ ^y = 0 . 


(24) 


The latter has solutions expressed in terms of Bessel functions. To see this we use the 
following known result. The equation 


x'^y" + (1 — 2p)xy' + 


2 2 2r I / 2 

q r X + [p — 


2 2 
n r 


y = 0 


(25) 
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has for real q the linear independent solutions 


yi = xPJniqx^), 


( 26 ) 


y2 = x^Yniqx''), 


(27) 


where Jn{x) and Yn{x) are the Bessel functions of the first and second type, respectively. 
We shall also use the following properties of the Bessel functions 


n 


X 


= 2 


"Tn— 1 ( 2 ^) 


n 


Jn{Y) '7n+l(2^)' 


X 


The same holds for Yn{x). Then, using (28) one gets 


Vi = yi 


p-nr ^_^Jn-i{qx^) 

+ qrx 


X 


Jniqx'^ 


(28) 

(29) 

(30) 


(31) 


1/2 = y2 


p-nr ^_;^Yn-i{qx'') 

+ qrx 


X 


On the other hand, using (30) one gets 

\p + nr 


yi = yi 


X 


— qrx 


Yn{qx'^ 


,_i Jn+i{qx^) 

Jniqx'^) 


(32) 


(33) 


1/2 = y2 


p + nr ^_^Yn+i{qx^) 

- qrx 


X 


Yniqx'^ 


(34) 


In order to hnd the Riccati solutions we can identify the parameters from comparison of 
(24) and (25) 


1 2 


ab 


3-6 


1 


r = 


n = 


(35) 


2 ’ ’ 3-iV r(2-i) ’ ' 2 ’ 3-5 

It is more convenient to work with equations (31) and (32) because p — nr = 0. Therefore 

r — \Jn—l(^qX ) 


( ^ y'l 1 

ui{x) = — = - 
ayi a 


qrx 


Jn{qx^ 


(36) 


^ 2 ( 2 :) = ^ = - 
ay 2 a 


qrx 


._iYn-i{qx'^) 

Yn{qx^) 


(37) 
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and thus 


ui{x) 


U2{x) 


ab 


T ( ^ / ab 

I-'’ r(2-5) 

-x~ ^ 


X 2 


r(2-<5) j^(2 
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3^V3-5V |r(2-5)| 


rX 2 


ab 


r(2 - 5) 


1-S ^|5| ( 3-5 \/r(2-5) 


a; 2 


r(2-5) 


afe 3-*^ ' 
ao ^ 2 


(38) 


(39) 


For positive b one gets an imaginary q parameter that turns the J and Y functions in the / 
and K Bessel functions, respectively. If we consider normal (nondivergent) initial conditions 
as a criterium for physical solutions than this selects the expressions containing the J and 
I functions. 


V. APPLICATION TO FRW BAROTROPIC COSMOLOGY 


Recently Faraoni [5] showed that the equations describing the FRW barotropic cosmolo¬ 
gies can be combined in a simple Riccati equation of constant coefficients. In addition, Rosu 
[6] discussed in some detail the cosmological Riccati solutions and used nonrelativistic su¬ 
persymmetry (generalized Darboux transformations) to get cosmological Riccati equations 
of nonconstant coefficients. Faraoni’s Riccati equation is 

^ + cH^ = -kc, (40) 


where H{r]) = jg the Hubble parameter (R is the scale factor of the universe) and 

T] is the conformal time, c is related to the adiabatic index of the cosmological fluid under 
consideration, c = |7 — 1. k = 0, —1,1 are the curvature indices of the FRW universes, 
plane, open, and closed, respectively. 

Applying the results of the previous section, for a = c, 6 = —kc, we get the following 
solutions for the S- modihed Hubble parameter: 

For k = 1 (the closed case), the q parameter is real and we get 







(41) 


(42) 


For k = —1 (the open case), the q parameter is imaginary and we get 

ff<->(.,;^) = 


r(2 - S) 


1-5 (3-5 V r(2-5) ^ 


3-5 . 
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2c 


^3^^3-5V r(2-5) 




(43) 
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(44) 




\/r(2-5)^ 


3-5 « 
2 




The case A; = 0 corresponds to 6 = 0, therefore it does not enter the present generalization 
in the sense that there is no change in the Riccati equation. 

The formulas (41-44) can be considered a generalization of the results obtained by 
Faraoni. Nondivergent initial data correspond to (41) and (43). Three-dimensional plots of 
these formulas are given in Figures (1) and (2). For 5 = 1 we get 5) = J- 1 / 2 / J+ 1/2 = 

cotan(c77) and H\~\rj]6) = /_i/ 2 //+i /2 = cotanh(c77) that correspond to the ordinary cal¬ 
culus. We mention that there are various works in the literature on the issue of geometric 
and physical interpretation of the fractional derivative and fractional integral, see, e.g., Pod- 
lubny [7]. In the cosmological case, the new parameter 5 is introduced in the cosmological 
evolution of the Hubble parameter as a consequence of applying a special fractional calculus 
to cosmological realms. In principle, as in statistical mechanics [8], the fractional calculus 
can be considered as the macroscopic manifestation of randomness. This has been argued 
to be so [8] when there is no dehnite time-scale separation between the macroscopic and the 
microscopic level of description and this could be the case of cosmology itself. 


H(eta;delta) for closed universe and c=1 



FIG. 1. Fractional Hubble parameter calculated according to the formula (41). 
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H(eta; delta) for open universe and c=1 



FIG. 2. Fractional Hubble parameter calculated according to the formula (43). 
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